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Abstract 

We study the effect of an electric charge in the middle of a ring of electrons in a magnetic 
field such as = 1/2. In the absence of the central charge, a residual current should 
appear due to an Aharanov-Bohm effect. As the charge varies, periodic currents should 
appear in the ring. We evaluate the amplitude of these currents, as well as their period 
as the central charge varies. The presence of these currents should be a direct signature 
of the existence of a statistical gauge field in the v = 1/2 quantum Hall effect. Numerical 
diagonalizations for a small number of electrons on the sphere are also carried out. The 
numerical results up to 9 electrons are qualitatively consistent with the mean field picture. 
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1 Introduction 



Recently, new experimental and theoretical developments have generated alternative viewpoints in 
the current understanding of the fractional quantum Hall effect. It seems the key idea is the notion 
of composite fermions, first introduced by Jain Q, which establishes a correspondence between 
the fractional quantum Hall effect at lowest Landau level filling v = p/[2mp ib 1) and the integer 
quantum Hall effect with p filled Landau levels. The mapping is achieved by attaching 2m flux quanta 
to each electron, which preserves fermionic statistics, and by first treating the statistical fluxes at 
the mean field level. This description has received considerable attention especially in the vicinity of 
v = 1/2. Experimentally, the system is characterized by a vanishing energy gap, and various physical 
properties strongly suggest the presence of a new Fermi liquid Q in this problem. On the theoretical 
side, the composite picture with m = 1 has been advocated by Halperin, Lee and Read as a powerful 
microscopic basis to understand these Fermi liquid-like properties in the vicinity of = 1/2 Q. 
The vanishing of the gap for the sequence of states aX v = p/{2p ±1) as p goes to infinity receives 
for instance a very simple interpretation in this framework since the composite fermions experience 
an average flux equal to 0o/p per particle, so the corresponding Landau level spacing goes as 1/p. 
If the statistical gauge field fluctuations are taken into account, some complications arise since a 
logarithmic divergence of the effective mass as p ^ +oo is predicted Q. There has been some 
trends in recent observation suggesting an effective mass enhancement but this question is not 
settled yet. This theory has also been tested by numerical diagonalizations on finite systems Q Q, 
and the composite fermion picture has been shown to predict for instance the right quantum numbers 
for the ground state and low- lying states , as well as the scaling of the ground state energy versus 
particle number Furthermore, a trial wave function inspired by these considerations provides a 
very good understanding of ground state correlations. We note however that this wave function is not 
merely the singular gauge transformation applied to the Slater determinant of composite fermions 
moving in zero external field, but it has to be improved by the combined effect of a short range 
Jastrow factor and a global projection on the lowest Landau level. 

In this paper, we suggest a different test for the composite fermion approach. The basic idea 
is connected to the fact that if an external electric field generates a spatially dependent electronic 
density, the average flux acting on the composite fermions is no longer vanishing everywhere. We 
wish to detect such a variation of the effective flux. In the experiment we propose, local density 
fluctuations are created in a ring oi v = 1/2 fermions, thanks to the presence of an electric charge 
localized in the middle of the ring. Due to the presence of this charge, the electrons have a trend 
to accumulate at the external edge of the ring if the central charge is negative, and to be repelled 
from the inner edge, which creates a non zero effective flux through the ring. A consequence of 
this non zero flux is the existence of observable persistent currents which vary periodically with the 
value of the charge in the center of the ring. At first glance, this phenomenon is reminiscent of the 
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Aharonov-Casher effect, whicli is expected to take place if a flux tube moves around a fixed charge. 
[^. This effect is a direct consequence of the Aharonov-Bohm effect and of Lorentz invariance, since 
the static electric field generated by the charge induces a non- vanishing magnetic field in the flux-tube 
co-moving frame. Experimentally, it has been observed with superconducting vortices in Josephson 
junction arrays However, we should emphasize here that the composite fermions do not carry 
the physical magnetic field which obeys Maxwell's equations, but rather a Chern-Simon flux, so the 
two situations are physically quite different. 

This paper is organized as follows. We first sum up the formalism to deal with the mean field 
approach of the ly = 1/2 quantum Hall effect. We then describe the geometry of the experiment, 
and propose a simplified geometry to get rid of the effects of curvature. In the absence of the central 
charge, one is able to solve exactly the eigenstates of the problem in the geometry of the experiment. 
In the presence of a central charge, the ring is electrically polarized. We propose a self consistent 
approach to find the electron density in the ring. Due to the presence of screening, the density 
fluctuations are localized in the vicinity of the edges of the ring. We finally obtain an order of 
magnitude for the average magnetic field through the ring. We calculate the current in the ring as a 
function of the central charge. We then present numerical computations of the effect on the sphere. 
After briefly writing the mean field theory on the sphere, we give numerical results for u = 1/2. A 
conclusion discusses the various results and mentions some open questions. 

2 Mean field theory at z/ = 1/2 
2.1 Gauge transformation at u = 1/2 

We first sum up the mean field theory treatment of the half field Landau level Q. It is possible to 
transform the original fermions in the magnetic field into new fermionic composite particles, with 
two flux quanta attached to them. The reason why attaching two flux quanta leads to fermions is 
that when one interchanges two particles with a flux (p attached to them, one gets a phase factor of 
exp i9 in the wave function, with 

9 = 7r{l + ^), (1) 

where 9 = ir refers to fermionic particles, (po = h/e is the flux quantum. The case (p = 2(j)Q thus 
corresponds to fermionic particles since in this case 9 is n modulo 27r. The flux tubes are fixed via 
a non local gauge transformation, implemented in the following way. The second quantized form of 
the kinetic term in the Hamiltonian reads 

^ = ^1 d'xV^+(x)(-i?iV + eA(x))2V'(x), (2) 

where M is the band electronic mass, e the absolute value of the electronic charge and V'^(x) the 
electronic field. The gauge transformation which realizes the passage from the electronic field ip'^^x) 
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to the composite fermionic field ip^{^) is given by 

^+(x) = V+(x) exp {-2i J (f^'arg{j^ - x')p(x')} (3) 

where 

;5(x) = ^+(x)V'(x) = V'^ (x)^c(x) = /5c(x), (4) 

and arg(x — x') is the angle between the vector x — x' and the x axis. The factor 2 in the exponential 
of stands for two flux tubes attached to the electrons to form the composite fermions. In terms 
of the composite fermions, the kinetic energy reads 

^ = 1^ / dxV^^ (x)(-iV + |a(x) - a(x))2^c(x), (5) 
where the statistical gauge field is 

a(x) = 2 / dx'^^^^i^pc(x'). (6) 



The interaction term is 

I r 1 



V = - J dxdx'i;(x - x') : p(x)p(x') ■ = 2 J d^dx.'v{-K - x') : pc'(x)/5c(x') :, (7) 
v{x) being the Coulomb interaction 



e2 



v{x) = --, (8) 

47reoer.|x| 



where is the relative dielectric constant of the material. For GaAs, = 12.6. 



2.2 Mean field theory 

In the mean field theory approach, the density of fermions p(x) is assumed to be constant. Since the 
average density uq is related to the field by the condition of half filling 

"» = 

the mean field value of the statistical flux exactly cancels the external magnetic flux, and gives a 
zero residual magnetic field. The Hamiltonian of the system of composite fermions is simply the 
Hamiltonian of a collection of free fermions, but with a renormalized effective mass M* 

Hm.f. = ^J ^+(x)(-iW)V(x)dx + V. (10) 

If B = WT, the value of the effective mass for GaAs is M* ~ 4M ~ 0.27 Me, where Me is the bare 
electronic mass. M* increases as the square root of the magnetic field for larger magnetic fields 
[^. If the density n(x) deviates from the average density no, then a residual effective magnetic field 
appears, equal to 

B - 20on(x) = 20o(no - n{x.)) = -20o<5n( x). (11) 
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3 Geometry of the experiment 



The electrons are confined on a small two-dimensional ring. The ring is the set of points x, such as 
ro < |x| < ro + L. The dimensions of the ring are of the order rg — 1/im and L ~ 0.1/xm. We look 
for the response of the system to an extra charge Q added in the center of the ring. As we shall see 
in the next section, one is able to solve exactly the Schrodinger equation of the electrons on the ring 
in the absence of the central charge and in the absence of Coulomb interaction between electrons. 
The wave functions are given in terms of Bessel functions. In order to simplify the treatment of the 
problem, we change the geometry, (see figure [l|). Instead of a ring, we shall use a rectangle of size 
L and R = 27rro. The x axis is chosen along the R side of the rectangle, and the y axis along the 
L side of the rectangle. We impose cyclic boundary conditions in the x direction, which means that 
ip{x + R,y) = ip{x,y) for the wave functions. The spectrum and the wave functions are much more 
simpler for the approximate geometry. Two types of potentials in the x direction will be investigated: 
an infinite square well and also a parabolic well, to understand the effects of a gradual fall-off of the 
electronic density as may be the case in experiments. 

4 u = 1/2 electrons on the ring in the absence of the central charge 

One is able to solve exactly the eigenvalue problem in the absence of the central charge, in the 
absence of a Coulomb interaction between electrons and in the ring geometry, as well as in the 
simplified geometry. The magnetic field is related to the density by the condition of half filling (|9|), 
so that the total effective magnetic field through the ring is zero, and the magnetic field in the hole 
inside the ring is uniform and given by 

B = 2<j)ono. (12) 
The number of electrons on the ring is given by 

N = -mL + 2r,)_ 

If we take B = 20T, the number of electrons on the ring is N = 1494. These N composite fermions 
feel a zero magnetic field, but feel the vector potential created by the flux tube in the hole of the ring. 



This is a typical Aharanov-Bohm situation [10| and one expects the presence of electronic currents 

ly 

(14) 



which are periodic in <p/(j)o, where if denotes the flux in the center of the ring. Its value is simply 

ip _ ■nr^B 



Similar situations have already been analyzed in different geometries |11] ||l^. The gauge is chosen 
such as A(x) = 74(|x|)e6), with 

Mr) = ^ (15) 
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if r > rg. One is left with a problem of free electrons in the vector potential A(r) with a Hamiltonian 

H = ^—{-inv + eAf. (16) 

Because of the rotational invariance of the problem, the wave functions can be chosen with a definite 
orbital momentum mh 

^l;{r,e) = e'^\{r). (17) 

The Schrodinger equation for x(r) reads 

^"(r) + ^ + (A:2 _ J_(^ + ^?)x{r) = 0, (18) 



which is a Bessel equation where we have set E = h?k'^ /2M* . The general solution of this equation 
is 

X(r) = AJ|„+^|(A;r) +i?l]^+^|(/cr). (19) 

The coefficients A and -B are determined by the normalisation of the wave function and by the fact 
that the wave function vanishes at the edges of the sample, so that xi^o) = xi^o + = 0- The wave 
vector k is found to be the solution of 

J|m+^|(/cro)l^„+^|(A;(ro + L)) - ■5^„+^|(fcro) J|^+^|(fc(ro + L)) = 0. (20) 

In order to obtain simpler equations for the energy levels and the wave functions, we give the form 
of the solution in the simplified geometry. We first treat the case (p = 0. In this case, the problem 
simply corresponds to free electrons on a rectangle. The wave functions are given by 



^mA^, y) = sin [n^y^ , (21) 

and the energy levels are 

The electrons belong to a Fermi sea, and the Fermi wave vector kp 'is approximately determinated 
by the relation 

(.3) 



where the number of fermions N is given by (|13D. Since L <^ R, the wave vector increment 2t:/R in 
the kx direction is much smaller than the increment vr/L in the ky direction, so that the Fermi sea 
can be viewed as a collection of channels labeled by the integer n > 0. The Fermi sea is drawn on 
figure 0. In the case B = 20T, the Fermi sea is made up of six channels. In the absence of a magnetic 
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flux (p through the hole of the ring, or if the magnetic flux is a multiple of the flux quantum (pQ, one 
can compute the electronic density of the Slater determinant 



lV'o>= n V'+JO), (24) 

m,n&F.S. 

where the fermionic quantum numbers m and n belong to the Fermi sea of figure ^ and |0) is the 
vacuum. The density profile is plotted on figure |3|. As expected, the density is zero on the edges 
and the density profile exhibits Friedel oscillations. The non uniformity of the density induces a non 
uniform electrostatic field, which modifies the one-particle states. Thus, due to finite size effects, (p^) 
is not the true ground state for the mean field approximation, whereas it would be the true ground 
state on an infinite plane. 

What happens if the flux ip in the hole of the ring is non zero? Let us first consider the case 
where (/7/(/>o is an integer. As we see from equation ([T8|), we can deduce the physics from the case 
(/5 = by replacing m by m + (/j/(/>o- The Fermi sea is translated in the reciprocal space by a factor 
A/cx = 2inp/R(j)Q along the kx direction. 

If ip/(t)Q is not an integer, the situation is somewhat different. One expects in this case the 



appearance of a current due to an Aharanov-Bohm effect |13]. We note v'/^o = + <^"^) with Am 
an integer and 5m a real number, such as \6m\ < 1/2. The energy E(m,n) reads 

27r vr 
E{m,n) = ^[{—)\m + Am + 6mf + (25) 

The Fermi sea is translated by the wave vector Ak^ = — 27rAm/i?, so that the wave functions become 



= ^_e^(— ^-)^-sinn^y (26) 

The Fermi sea becomes unstable as \m\ = 1/2. To see this, consider 2A''o + 1 fermions in a given 
channel n. If \6m\ is less than 1/2, the fermions have the possibility to have their orbital quantum 
numbers in the interval [—Am — A'^o, — Am + Nq], or in the interval [—Am — A'o ± 1, — Am + Nq it 1]. 
The condition for the latter configuration to be stable is that it has a lower energy than the former, 
namely 

— — (— )^ V (m + Am + 5mf < ——(—f V (m + Am + 5m)^, (27) 

m=—Am—NQ m=—/\m—No±l 

which is satisfied for |5m| < 1/2, independently on Nq. The Fermi sea of 2Ao + 1 fermions in the 
channel n is thus in the interval m E [Am — NQ,Am + Nq] if —1/2 < 6m < 1/2. What is the 
Aharonov-Bohm current in the ring for a non integer value of ip/(j)Q ? The quantum mechanical, 
gauge invariant current operator reads 

j W = ^(V^W - V'VV^) + ^A|^|2. (28) 
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In the simplified geometry, the current has only a component parallel to the x axis. The contribution 
of the first term in (|2^) is 

— — m-Am, (29) 
and the contribution of the term proportional to the vector potential reads 

(Am + (5m). (30) 



M* R 

The total current is thus 

jx{m,n) = — — {m + 6m)\'>pm,n\'^- (31) 
After a summation over the Fermi sea, we obtain the total current 

Jx{y)= J2 j^i^^^) = ^-^^^ IV'm,n|^ (32) 

{m,n)eF.S. {m,n)£F.S. 

If 5m = 0, the flux in the hole of the ring is a multiple of the flux quantum, and everything happens 
as if the fermions would not see the flux. If 5m 7^ 0, a current exists. The intensity of the current is 
given by 

I = -e I Jx{y)dy = - ^^^^ 5m. (33) 







The maximum value of the current is 



^max - 2M*R^- ^^^^ 
The numerical value of Imax is 5.8nA for B = 20T. 

5 u = 1/2 electrons on a ring in the presence of the external charge 

In the presence of a positive central charge, the electrons are expected to accumulate near the internal 
edge of the ring and a depletion of electrons is expected at the external edge, leading to a charge 
transfer from the external to the internal edge. This charge transfer induces a deviation from the 
V = 1/2 value and, if ones applies the ideas of the mean field theory, the total effective flux, due to 
the external magnetic field plus the statistical gauge field is no longer zero, so that an effective flux 
penetrates through the sample, creating an Aharanov-Bohm current which is periodic as a function 
of the central charge. 

The first step is to evaluate the screening of the central charge by the electron gas on the ring. 
Many levels of approximation are possible. The most accurate approximation is to take into account 
that the one particle Slater determinant (p^) is not the true ground state of the fermions, due to 
Friedel oscillations induced by the edges. This can be implemented in a recursive way by starting 
from the state (|2^), computing the electronic density, deducing from it the electrostatic field, and 
reiterating, that is to compute the corrections to the one-particle states in the presence of the Priedel 
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oscillations, and re-compute the electronic density. The effect of the charge is then treated at the 
linear order, with the full response function. Within these approximations, one expects the density 
of electrons to increase near the internal edge of the ring, and to decrease near the external edge. 
Moreover, this average density profile should be modulated by Friedel oscillations as well as oscilla- 
tions at the Thomas-Fermi wave vector. In our case, as we shall see, the Thomas-Fermi wave vector 
is greater than the Fermi wave vector. We shall not use this refined approximation scheme because 
of its computational complexity. Since we focus only on the charge transfer from one edge of the 
ring to the other, and not on the details of the variations of the density profile, we do not take 



into account the existence of Friedel oscillations. We take the state (24) as an approximate ground 
state and we look for the effects of the screening of the central charge. In other words, we hope the 
density-density response functions are not very different if we replace the true Hartree-Fock state by 
the approximate one given by (p^). 

At this stage, two approximations are possible. The most refined one consists in taking into 
account the full response function of the fermions. This shall be done in section 5.4. A more 
approximate treatment is to treat the screening in the Thomas-Fermi approximation, which is the 
aim of section 5.3. As we shall see, the two approximations lead to similar results as far as averaged 
quantities are concerned, namely the average effective flux penetrating through the ring due to the 
presence of the central charge. 

5.1 Screening in two dimensions 

As an introduction to the problem of screening in a finite geometry, we treat the case of the screening 
of a single charge in an infinite two-dimensional gas of electrons. The solution of this problem will lead 
to the expression of the Thomas-Fermi wave vector qxp in two dimensions. The case of the screening 



of a charge in three dimensions is treated in reference |14] with the Thomas- Fermi approximation. 



We do nothing but transpose the argument of |14] to the case of an electron gas constrained on a 



two-dimensional layer, with three dimensional Coulomb interactions. The potential created by the 
external point charge located at the origin is noted c/)^^*. We note /O*""^ the variation of density induced 
by the presence of the extra charge in the gas of electrons, plus the uniform background of positive 
charges, (j) is the total potential created by the extra charge, the electrons and the background of 
positive charges. Since the two-dimensional Fourier transform of l/|x| is 27r/|q|, we have 

Jnd( \ 

<A(q)-0-*(q) = ^^. (35) 

The dielectric constant is defined by (/)'^^'*(q) = e(q)i?i)(q), and we assume that p*"'^(q) = x{^)4'{^}■ 
We thus obtain 

e(q) = 1 - (36) 
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If the total potential is a slowly varying function of the position, one can define 

^(k) = 1^ - e</.(x), (37) 

so that the local distribution function reads 

r dk 1 

~ J (27r)2 1 + exp {/3(h2k2/2M* - e0(x) - fj,)} ' ^ ' 

and the density of the background positive charge is 

"""^^^ = / (27r)2 1 + exp {/3(n2k72M* - fi)} ' ^^^^ 

Thus, we can write the induced density of electrons as /9*"'^(x) = — e{n''(/x + e(/)(x)) — n^{ij,)}. If ^ is 
small, p*"'^(x) = — e^5n°/3/i(/)(x), so that 



The dielectric constant reads e(q) = 1 + ^tf/IqI, where qte is the Thomas-Fermi wave vector 

h"^ 47reoe 

If one adds an extra charge Q at the origin, 



9TF = T2-7rr^- (41) 



A Fourier transform yields 



with 



(42, 

F{z) = Mu)^^du. (44) 

The interaction is screened if x 2Tr/qTF- The numerical value of the Thomas- Fermi length Xtf = 
217 /qtf is Atf = 15.6nm. 

5.2 Linear response approach 

The aim of the section is to present the real space linear response formalism. From a numerical 
point of view, this means that wc must discretizc the radial coordinate at a scale smaller than the 
Thomas- Fermi length Xtf- Within the linear response, the density variation is linearly related to 
the local potential, but in a non local way: 

<5p(x) = x(°^ (x, x')Vioc{x')dx', (45) 

JO 



9 



where x E [0, L] is the radial coordinate , with the origin taken at the interior edge of the ring. Notice 
that the conservation of charge carriers imphes that 

i-L 

X^°^(x,x')(ix = 0. (46) 



The local potential is the sum of the potential created by the electrostatic field of the central charge 
Vextix), plus the potential Vind{x) induced by the electrons on the ring. Since the distance L between 
the two edges is small (L = 0.1/im), we linearize the Coulomb potential created by the central charge 
in the vicinity of the ring, which leads to the following expression of Vextix)- 

Vextix) = ^ --^{x - -) (47) 

On the other hand, the induced potential has the expression 

„2 



Vindix) 



■ /^(l + —)6p{x')dx' (48) 
Jo ro 

(49) 

(1 + fi? + (1 + 'i? - 2(1 + + fj) cosO 

Using the notations u = I + x/vq and v = 1 + x' / rg , the angular integral takes the form of an elliptic 
integral 

dO . r/"^ dO 4 



'-TT V-u^ + t;2 _ 2uv COS 9 Jo \/{{u + v)'^ - 4uv cos^ 6 u + v \u + v 
Since u and v are close to unity, the elliptic integral can be approximated as 



A-l^ ■ (50) 



K 

\u + V I \ \u — v\ 

which leads to 

Vindix) ^ r In (y^^ ) 5p{x')dx' (52) 

2TTeoer Jo \\x - x" ' 

Notice that the factor 8ro does not come into account since 

rL 



6p{x)dx = 0. (53) 

JO 

The equation ([45D reads 

Sp = X^^\Vext + Vind5p), (54) 

where the quantities are understood as matrices for x'^ ^'iid Vmd and vectors for 5p and Vio^- A 
discretization of the radial coordinate has been assumed, and we used a parabolic approximation for 
the integrals. A special attention has to be paid to the logarithmic divergence of ( |5^ ) which has to be 
integrated explicitly using a parabolic approximation. Equation (|5^ ) can be inverted into the form 
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which has the well known form of a R.P.A. resummation. Numerically, we use a Gauss-Jordan 
method for inverting the linear system (|52|). The size of the matrix to be inverted is about 1000, 
which corresponds to 1000 points for the discretization of the interval [0, L]. Since L ~ 0.1/im, the 
condition L/1000 ^ Atf is well respected. We now examine successively two levels of approximation 
for the response function x^^\x-,x')- 

5.3 Thomas-Fermi approach to the response function 

In the Thomas-Fermi approach, the response function is purely local: 

xi?^(x,x') = --^(5(x-x')--) (56) 

Using this form of the response function, we inverted the system (54). The resulting density profile 
is given on figure In terms of the effective flux through the sample (flux of the magnetic field plus 
flux of the statistical gauge field) , one is interested in 

(t){x) = -A^(t)o{ro + L/2) r 6p{x)dx. (57) 

Jo 

The function (j){x) is plotted on figure I As we shall see in section 5.5, the total current is given in 
terms of the mean flux through the ring 

- 1 



L 







(P{x)dx (58) 



We find that is such that a charge Q = 14e~ is required to induce one flux quantum on average 
through the ring, which means that the periodicity of the currents as a function of the central charge 
is 14 electrons within the infinite square well potential and the Thomas Fermi approximation. 

5.4 Screening of the central charge with the full response function 

We calculate the full response function using the first order perturbation theory in the potential 
induced by the central charge. We also use the simplified geometry. The matrix elements of the local 
potential V/oc on the basis of function (^) are 

Vn%' = J^J^ sin (^^rax^ Viocix) sin (^j^n'x^ dx, (59) 

so that 

^oc = E E Kn'i^t^,ni^m,n' ■ (60) 



The first order corrections to the state \ipo) of equation (l24| ) are 



I V'O) + E E /^"'"f^ N V'm,nV'm,ri' I V'o) • (61) 
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The second-quantized form of the electron density operator is 

p(^^y) = E Ee"^^'"'~'"^^''si" il^y) (^^'y)<nV'n^',n'iV'o>. (62) 

m,m' n,n' 

One can readily calculate the average of p{x,y) in the presence of the central charge, which leads to 

"m n,n' \ ' / \ j / 



(1 — vP{m, n)) sin ^^^^^ ^™ (l^'^) ' ^^^^ 



where nP^n, m) is unity if (m, n) belongs to the Fermi sea and zero otherwise. The summation over m 
can be evaluated by noticing that the Fermi sea is made up of a collection of nested one dimensional 
Fermi seas (one for each channel), and we obtain 

5(^(^,2/)) = ^ { E E ^^^2 : ^"2' " {l-y) {l-'y) 

n=l n—n+1 

N +00 / 9 9 / 

n=l n'=Af+l 



^ / 2 2 / \ / \ 

+E E ^^-HH'nPT"-'^- ^''^ 



In this expression, v is related to the Fermi energy by 

2M*L2 ^ 
and is the integer such as G [A^, + 1[. We deduce from ( |64|) the response function x^'^H?/) ?/')• 

X^^ = ^ E E ^^37^7^ (2/-^) + E E T^^37^Al,„, (y,y) 

rt=l rt'=rt+l n=l n'=Af+l 

(66) 



with 

The symbol IFW stands for "infinite square well". With the form ( |66| ) for x^^\ it is straighforward to 
compute the density of electrons, using the matrix relation (^). The profile of the electron density 
variations induced by the central charge is plotted on figure ^ We notice the existence of oscillations 
at the Thomas-Fermi wave vector, which were absent in the calculation using the Thomas-Fermi 
approach to the response function. However, these oscillations do not affect too much the averaged 
quantities under interest. The total fiux through the sample 0(x) is plotted on figure and has the 
same shape as the flux computed in the Thomas-Fermi approximation. The mean flux through the 
ring is such as 18.2 e~ are required to produce one flux quantum through the ring within the infinite 
square well approximation and within a calculation including the full response function. 
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5.5 Electrons in a parabolic potential well 

In experiments, the edges of the ring may not be so weh-defined, and the electronic density will 
exhibit a gradual fall off to zero. This section is devoted to model the effects of ill-defined edges. To 
do so, we use the same techniques as before, and use a parabolic confining well instead of an infinitely 
deep well. Within the mean field treatment, the Hamiltonian of the quasiparticles reads 

rP' 1 

H = + -M*uj^y^. (68) 
2M* 2 ^ ^ 

As before, we start from a Slater determinant (24), with one electron wave functions 



where the inverse length scale /? is 




ip^,nix,y) = ^lt]' -^exp ( W„(/3y)e-™1r- (69) 



, 2n , 



^-r-^ (TO, 

and Hfi are the Hermite polynomials. The radial density of the Slater determinant is plotted on figure 
where the strength of the harmonic potential has been adjusted such as 70% of the electronic density 
lies within the one micrometer wide ring. Friedel oscillations are visible. From a computational 
point of view, the Hermite polynomials may reach huge values. It is thus more practical to calculate 



recursively Hn = Hn/V^^^nl, with the following recursion relations 



Hn+liu) = J —uHn{u) - J—^Hn-l{u), (71) 



and Hq{u) = 1, Hi{u) = ^/2u. The same approach as in section 5^ can be carried out and the 
response function is 

\n=Qn'=n+l n=On'=N+l 

(72) 

with 

K,n'{y,y') = Hn{f3y)Hn'mHn{f3y')Hn'{Py'), (73) 



and v is defined by Ep = (i^ + l/2)huj. Using similar techniques as in section bA, we get the radial 



density variations (figure P) as well as the induced flux (figure 10). We come to the conclusion that a 
central charge approximately equal to 9 electrons is required to generate one flux quantum through 
the ring, which suggests that the Fermi sea in a parabolic well is more sensitive to the central charge 
as in the infinite square well. This is due to the fact that the fraction of the electronic density outside 
the O.lfim ring experience is closer to the central charge (as far as the inner edge is concerned) and 
thus experiences a stronger potential. We conclude to the importance of whether there is a sharp 
edge or a graduate fall-off to zero of the electronic density. In the former case, one flux quantum 
could be generated through the ring with a smaller charge than in the latter case. 
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5.6 Currents induced by the charge 

The magnetic flux through the ring generates permanents currents, in the same way as a flux through 
the hole of the ring generates currents. Let us call 6 A the variation in the vector potential due to 
the presence of the charge. The variation in the Hamiltonian due to the shift in the vector potential 

6H = ^(2A.<5A + 6 A') + ^^A.V. (74) 
Applying first order perturbation theory, we obtain the variation of a given energy level 

{SH) = ^(A.5A) - ^m(^). (75) 

^ ' M* M* r ^ 

We now make the approximation that the spatial dependence of 6 A and 5A{r)/r are averaged on 
the ring, that is to say 

6A{r) ^ 



(76) 



where (p is the average magnetic flux induced by the central electric charge. The flrst term in 6H is 
simply an energy variation independent on n and m. This term is then dropped. Thus we get 



eh 

-m — - 

so that the total energy reads 



{SH) = -— m— ^, (77) 



E{m,n) = ^-^{m-^f-^-^{^f. (78) 

If we apply the result of section 4 about Aharonov Bohm currents, we find that the currents are 
periodic in (p. The maximum value of the current is given by (p3) and the period is such as cp = (po- 



6 Numerical approach in the spherical geometry 

We now make use of the spherical geometry in order to perform numerical computations with a small 
number of electrons on a sphere of radius R. A magnetic monopole is put at the center of the sphere, 
creating a total magnetic fiux through the sphere equal to 2S(pQ. The Dirac's monopole quantization 
requires 25 to be an integer. An electric charge Q is put at the north pole of the sphere. The extra 
charge is treated as a classical, point charge. We first generalize the mean field theory argument to 
the spherical geometry case. In a second step, we present numerical computations for a small number 
of electrons on the sphere for various filling fractions. 
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6.1 Mean field theory at u = 1/2 

We adapt the argument of the mean field theory, previously established in the geometry of the ring, 
to the spherical geometry. The filling fraction 

N -I 



■ 2S 

is chosen to be 1/2 in this section. The non local gauge transformation leads to a statistical field 
which generates a flux equal to —2(j)Q{N — l) because there is no statistical flux coming from a particle 
onto itself. The mean field Hamiltonian simply corresponds to fermions on the sphere in the absence 
of a magnetic field. If one neglects the Coulomb interactions, the Hamiltonian is diagonal on the 
basis of the spherical harmonics Yi^m{0,ip) normalized such as 

Yi^U9,^)dn = l (80) 

The wave function is simply 'i/'i,m(^) ^) = ^,m(^) f)/^- The single particle states are labelled by the 
positive integer / and the integer m such as — / < m < L The energy of a state labelled by {l,m) 
reads 

E(l, m) = ^1(1 + 1). (81) 

In the presence of a negative (positive) electric charge at the north pole, a depletion (accumulation) of 
electrons arises at the north pole, and an accumulation (depletion) arises at the south pole. The mean 
field model in the presence of the charge at the north pole consists of electrons in a zero magnetic 
field plus a flux tube (j) penetrating through the south pole and emerging at the north pole. The 
problem of quantum motion around a flux tube (p = a(/>o in the spherical geometry was considered in 
[ p!6[ . We now rederive the solution. This problem is non perturbative in 0, since inflnite quantities 
appear in the first order perturbation theory in (j). This is why we adopt an algebraic approach. The 
flux tube (j) is absorbed in a gauge transformation, leading to multivalued wave functions 

ipir, e + 2TT) = e^'^^^{r, 6). (82) 

The eigenvalues of Iz are thus quantized hy Iz = m + 27r7, m being an integer. The eigenstates of 
the Hamiltonian are also eigenstates of P, since 

A priori, it is not obvious to produce a basis of commun eigenvectors to P and Iz-, since the rotational 
invariance is broken by the presence of the flux tube. Nonetheless, as we shall see, it is possible to 
diagonalize simultaneously P and Iz- The spherical representation of the kinetic momentum algebra 
is Q: 

= =-(|+"o'«^) m 
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r = e-*'^(-^ + icot^^ ) (85) 



de dip, 

The operator algebra gives rise to two ladders of states. The (— ) ladder corresponds to states 
descending, by the repeated action of Z~, from the highest weight state such as l~^ipQ~^ = 0. 

The states tpQ have the form 

^i-\9,ip) = f<-\e)e'^'+^^^, (87) 

with I an integer, and is such as 

+ ' "^^^ e'^'^^^'"/^-' (0) = 0, (88) 

so that /^"H^) « (sin^)'+^, with Z + 7 > 0. We take 7 g]0, 1[, and I = 0,1, ... For a given I, we can 
produce I + 1 descending states by the repeated action of /~. For these states, Iz = 7, 1 + 7, .... Z + 7 
and F = (Z + 7)(Z + 7 + 1). The (+) ladder correspond to states ascending, by the repeated action 
of from the highest weight state i/j^^^ such as I'i^^^^ = 0. if)^^ has the form 

i^t\0M = f^^\0)e'^'^^^^, (89) 



with / an integer, and is such as 



+ icot^^ ) e^('+T)^/W(^) = 0, (90) 



80 dip, 

so that f^^\e) oc (sin^) with / + 7 < 0, so that |^| < —1. The repeated action of l'^ produces 

\l\ states with = ^ + 7, / + 7 + 1, -1 + 7 and 1^ = + 7)(/ - 1 + 7). 

We now enumerate the first states. With I = E = 7(7 + 1) and the degeneracy g is 1. With 
Z = -1, ^ = (1 - 7)(2 - 7) and g = 1. With I = I, E = {I + 7)(2 + 7) and g = 2. If / = -2, 
E = (2-7)(3-7) and g = 2. Iil = 2,E = (2 + 7)(3 + 7) and g = 3. If / = -p, E = (p-7)(p+l-7) 
and g = p and ii I = p, E = (p + 'y){p + 1 + 7) and g = p + 1. 

We now turn to the calculation of permanent currents, since this is the quantity we shall compute 
using numerical diagonalizations for a small number of electrons. We distinguish between two cases: 
pip-l) < N < p^ and < N < p{p + 1). First, if pip - I) < N < p^, we note X = N - p{p - 1). 
The total energy is given by 
p-i 

Etot = ^{{m-l + |7|)(m + I7I) + (m - \j\)im + 1 - I7I)) m + (p - 1 + \^\){p + (91) 

m=l 

= \p\p-lf+p{p-l)\ + i{2p-l)\-p{p-l))\j\ + {\ + p{p-l))j^ (92) 
for 7 G [-1/2, 1/2]. If p2 < AT < p{p + 1), we note \' = N - p^ and, for 7 G [-1/2, 1/2], we obtain 

Etot = \p\p - I? +p\p - 1) +P{P + 1)A' + - {2p + 1)A') |7| + if + y) 7^. (93) 
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The current dEtot/dj is discontinuous for half integer values of 7, as for the ring, but further discon- 
tinuities appear for integer values of 7. The jump in the current at 7 = is 

^(0+) -^{0-) = 2 {{2p - 1)A - p{p - 1)) , (94) 

or 

^(0+) - ^(0-) = 2 (p2 - (2p + l)y) . (95) 

a7 a7 ^ ^ 

The first term corresponds to the case p{p — 1) < N < p"^, and the second case to < < p{p + 1). 
The jump of the current at 7 = 1/2 is 

dEtot 1+ dEtot V 
d-f 4 ' d-i 4 
or 

^^2 )-^(2 ) = -4p(P-A). (97) 
6.2 Numerical procedure 

We assume that the magnetic field is large enough to neglect the excitations from the first Landau 
level to higher Landau levels. We propose to diagonalize numerically the Coulomb interaction inside 
the Hilbert space generated by the set of the occupations of the lowest Landau level. The quantum 



(- ) = -4pA, (96) 



Hall effect on the sphere has been studied by Haldane [15|. He proposes a set of coherent states 
which span the lowest Landau level Hilbert space. From this set of coherent states, we can extract 
the following basis of the lowest Landau level Hilbert space 

= (x|a) = iV„(sin^)"(cos^)2^-e^'^(^-°), (98) 
where is chosen in such a way that 

y'|0«(x)|2dx = l. (99) 

The label a indexes the orbitals of the lowest Landau level and runs from to 25. One needs to 
compute the matrix elements of the Coulomb interaction. The interaction Vi between the electrons 
on the sphere and the classical charge located at the north pole is a one-body operator and the 
interactions V2 between the electrons on the sphere is represented by a two-body operator. Using the 
second quantization, it is straightforward to derive expressions for the matrix elements of Vi and V2 
on the basis of the states with all the possible occupation numbers of the lowest Landau level. We 
do not give here the details of the calculations. We used the the results of |2C] relating the matrix 



elements of the Coulomb potential to Clebsh-Gordan coefficients. One is left with a symmetric 
matrix to be diagonalized using a Lanczos method in each sector of the z component of the angular 
momentum (since the Hamiltonian is invariant under any rotation around the z axis). The ground 
state is used to calculate the expectation value of the one-body current operator. The total intensity 
is obtained after an integration of the current j{9) over the angular coordinate 9. 
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6.3 Results 

We diagonalized the Coulomb interaction at half filling with N = 8 electrons (figure ^) and A'^ = 9 
electrons (figure |l3|). Obviously, the intensity is not periodic as the central charge varies. This is due 
to the fact that only a small number of electrons is present in our numerical diagonalizations. As 
the central charge is tuned, level crossing occurs, namely the z component of the angular momentum 
of the ground state is varied as the central charge is varied. Notice that these level crossings are 
allowed because they correspond to crossing of levels in different sectors of the Hilbert space. Due 
to the small number of electrons, as the central classical charge is increased, one reaches a point 
where no more crossing occurs since the Hilbert space, restricted to the lowest Landau level is finite. 
Above this point, the intensity saturates. If the central charge is large and positive, the electrons are 
attracted at the north pole and the current is thus positive. By contrast, if the classical charge is 
large and negative, the electrons are repelled from the north pole, and mainly located on the south 



hemisphere, and the current is negative (see figure 14). However, if the classical charge is not too 
important, we observe discontinuous variations of the intensity as a function of the central charge, 
which is very much reminiscent of the previously calculated mean field behavior even though we were 
not able to related the position of the discontinuities to the ones of the mean field model. We checked 
that the discontinuities in the current correspond to level crossings at the bottom of the spectrum 
in the sense that each jump in the current is associated with the fact that the z component of the 
angular momentum of the ground state changes. We conclude that the mean field scnario for the 
intensity as a function of the central charge is qualitatively supported by our numerical calculations. 
At the quantitative level, it seems rather difficult to conclude since we can only carry out numerical 
diagonalizations for a small number of electrons. In particular, we cannot conclude on the existence 
of periodic variations of the intensity as a function of the classical charge since the available sizes are 
too small. 



7 Conclusion 

We have shown that a charge in the middle of ring of = 1/2 electrons induces currents which vary 
periodically as a function of the central charge. The mechanism for the generation of these periodic 
currents involves a polarization of the electron liquid on the ring. If the extra charge is positive, 
a positive density variation appears on the interior edge of the ring, and a depletion of negative 
charges appears on the external edge of the sample. These charges screen the field of the central 
charge. The appearance of currents on the ring is mediated by the variation of the Chern-Simons 
gauge field at mean field level. The presence of polarization charges on the edges on the sample 
induces a non-zero average fiux through the ring. By increasing the charge, one should be able to 
produce periodic currents in the ring. The order of magnitude of these Aharonov-Bohm currents is 
such that they can be measured [^]. The amplitude of the current fiuctuations is typically 2.0 nA 
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for a ring of radius 1 and transverse dimension 0.1 fim and a magnetic field of 20 T. If one can 
produce a continuous charge in the middle of the ring, one should be able to find experimentally the 
period of the phenomenon. Our mean field calculations show that the period is of the order of 14 
electrons (calculated with the full response function) with an infinite square potential. However, the 
periodicity of the phenomenon is sensitive to the existence of a graduate density fall-off at the edges, 
as is expected to occur in experiments. In the case of a parabolic well with 30% of the electronic 
density outside the 0.1 micrometer ring, we found a 9 electrons periodicity. Residual currents should 
also appear in the absence of the charge, because of the presence of a magnetic flux through the ring. 
The experimental observation of such currents should be a direct test of the existence of a statistical 
gauge field in the v = 1/2 quantum Hall effect. We also carried out numerical computations in the 
spherical geometry. The mean field picture at = 1/2 predicts discontinuities in the intensity plotted 
as a function of the extra charge and also the existence of periodic variations. In our simulations, the 
current is not periodic as a function of the charge at the north pole, because of a too small number of 
electrons. However, we find the existence of discontinuities, reminiscent of the predictions of the mean 
field treatment. We conclude to the qualitative consistency between our numerical diagonalizations 
and the mean field theory. 

Our analysis shows that the type of experiment we consider would be a good test for the mean 
field theory of the = 1/2 quantum Hall effect since the presence of an extra charge allows to explore 
excited states properties, and this in a way which goes beyond linear response theory, since it involves 
level crossings and non-trivial reshuffling of the many body ground state as the external charge is 
varied. 

Finally, the case = 1/3 would also be an interesting case to examine and the physics might 
be quite different from the ly = 1/2 case investigated here. It would be especially interesting to 



understand the interplay between uncompressible bulk excitations and compressible edge states ^1 
as the central charge is varied. This question will be addressed in a future work. 

The authors acknowledge L. Levy for stimulating discussions which encouraged us to carry on 
this investigation. R.M. acknowledges J.C. Angles d'Auriac and F.V. de Abreu for assistance with 
the numerical work. 
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Figure captions 

Figure 1: 

(a): annulus geometry, (b): simplified geometry. 



Figure 2: 

Fermi sea of quasiparticles at the mean field level in the simplified geometry, kx is a multiple of 2tt / R 
and ky is a multiple of tt/L. 



Figure 3: 

Electronic density of the Slater determinant ( p^ ) p{x) = |(x|'(/'o)P- The electronic density vanishes 
at the edges and exhibits Friedel oscillations. The fermi sea contains 6 channels and 1494 electrons. 
The density is plotted in fim~'^ and the radial coordinate in fim. The electrons are confined in an 
infinite square well. 



Figure 4: 

Profile of the electronic density variations induced by the central charge in the Thomas-Fermi ap- 
proach to the response function. The computation has been made in the simplified geometry. The 
central charge is Q = —le~. The density is plotted in fim"'^ and the radial coordinate in fim. The 
electrons are confined in an infinite square well. 



Figure 5: 

Flux through the ring in the Thomas-Fermi approach to the response function. The central charge 
is Q = —le~. The flux is plotted in units of the flux quantum (pQ and the radial coordinate in fim. 
The electrons are conflned in an inflnite square well. 



Figure 6: 

Proflle of the electronic density variations induced by the central charge with the full response 
function. The computation has been made in the simplifled geometry. The central charge is Q = 
— le~. The density is plotted in ij.m"^ and the radial coordinate in /xm. The electrons are conflned 
in an infinite square well. 



Figure 7: 

Flux through the ring with the full response function. The central charge is Q = — le~. The flux is 
plotted in units of the flux quantum (j)Q and the radial coordinate in ^m. The origin of the radial 
coordinate is chosen at the inner edge. The unit length is the micrometer and the density is measured 
in micrometers"^. The electrons are conflned in an inflnite square well. 
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Figure 8: 

Radial density of the Slater determinant with a parabolic potential well. The Fermi sea contains 8 
channels and the strength of the parabolic potential has been adjusted such as 70% of the density 
is located on the 1 micrometer wide ring. Friedel oscillations are visible. The radial coordinate is 
counted in micrometers from the inner edge. The density is measured in micrometers"^ . 

Figure 9: 

Variations of the electronic density induced by a central charge equal to minus one electron, with 
a parabolic confining potential. The origin of distances is the center of the ring (minimum of the 
potential). The distances are measured in micrometers and the density in micrometers"^. 

Figure 10: 

Flux through the ring induced by a central charge equals to minus one electron with a parabolic 
confining well. The distances are measured in micrometers and the densities are in micrometers"^. 
The origin of the radial coordinate is the center of the well. 

Figure 11: 

Variations of the electronic current as a function of the central charge in the mean field approach. 
The period of the variations depends of the type of approximation (Thomas Fermi or full response 
function) and on the nature of the confining well (see the text). The residual current at Q = 
originates from the Aharanov-Bohm currents induced by the presence of the flux in the hole of the 
ring. 

Figure 12: 

Current as a function of the classical charge located at the north pole. The intensity is / = / j{9)d9. 
The radius R of the sphere is taken to be 1. The charge — Q is plotted in units of the electronic 
charge e. The current is plotted in units of 2eh/M*. The system contains 8 electrons. 

Figure 13: 

Current as a function of the classical charge located at the north pole. The intensity is / = / j{9)d9. 
The radius R of the sphere is taken to be 1. The charge — Q is plotted in units of the electronic 
charge e. The current is plotted in units of 2eh/M*. The system contains 9 electrons. 

FigTuc 14: 

Current circulation on the sphere for a large positive classical charge (a) and a large negative classical 
charge (b). In (a) the electrons gather at the north pole, due to the attractive Coulomb interaction 
between the electron gas and the classical charge. The current in this case is clockwise. In (b), the 
electrons are repelled from the north pole and the circulation of the current is counterclockwise. 



23 



3500 I 1 1 1 1 1 1 1 r 




Central charge = -1 electron 

0.5 I 1 1 1 1 1 1 




5 I I I I I I I I I I I 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 

Distance counted from the inner edge of the ring (micrometers) 



Induced flux variations witli Q = - 1 electron 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 

Distance from the inner edge (micrometers) 



Radial electronic density of the Slater determinant with a harmonic confining well 




-0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

radial coordinate 



Variations of the radial electronic density induced by a central charge = - 1 electron 




-0.15 -0.1 -0.05 0.05 0.1 0.15 

radial coordinate (micrometers) 



Flux thourgh the ring with a parabolic well 




-0.15 -0.1 -0.05 0.05 0.1 0.15 

radial coordinate (micrometers) 



8 



9 electrons on the sphere at half filling 



<x>oooooocoooooooooooooocooooooooooc<xxx>^^ 

6 - 



4 - 



2 - 




-8 I \ \ \ \ \ 1 

-6 -4 -2 2 4 6 

-Q 



